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jy^ ■ ABSTRACT 

I Context. The Red MSX Source (RMS) survey is an ongoing effort to return a large, well-selected sample of massive young stellar 

objects (MYSOs) within our Galaxy. 2000 candidates have been colour-selected from the Mid-course Space Experiment (MSX) point 
source catalogue (PSC). A series of ground-based follow-up observations are being undertaken in order to remove contaminant objects 
(ultra-compact HII (UCHII) regions, planetary nebulae (FN), evolved stars), and to begin characterising these MYSOs. 
Aims. As a part of these follow-up observations, high resolution (~1") mid-IR imaging aids the identification of contaminant objects 
which are resolved (UCHII regions, PN) as opposed to those which are unresolved (YSOs, evolved stars) as well as identifying YSOs 
Q ■ near UCHII regions and other multiple sources. 

V-^ ' Methods. We present I0.4/jm imaging observations for 346 candidate MYSOs in the RMS survey in the Southern Hemisphere, 

I primarily outside the region covered by the GLIMPSE Spitzer Legacy Survey. These were obtained using TIMMI2 on the ESO 

, 3.6m telescope in La Silla, Chile. Our photometric accuracy is of order O.OSJy, and our astrometric accuracy is 0.8", which is an 

improvement over the nominal 2" accuracy of the MSX PSC. 

Results. Point sources are detected in 64% of our observations, which are expected to be either YSOs or evolved stars. 24% contain 
only sources of extended emission, which are likely to be either UCHII regions or, in a few cases, PN. This is confirmed by comparison 
with radio continuum observations. We find that, as expected for a dusty HII region, the strength of 10.4yum and radio continuum 
emission is related. The remaining targets (12%) result in non-detections. While for 63% of our targets we detect only one mid- 
infrared source, 25% show multiple sources. In these cases, our observations will allow the apportioning of the flux from larger beam 
measurements between the different sources. Within these multiple source targets, we find some point sources on or near UCHII 
regions. Our improved astrometric information will allow more accurate targeting of spectroscopy, which will be used to identify 

■ unresolved sources in cases where it is not clear whether they are YSOs or evolved stars. 

■ Key words. Stars: Formation - Stars: Pre-Main Sequence - HII Regions - Planetary Nebulae 

o 

. ^ 1 . Introduction star, but has not yet begun to ionise the surrounding molecular 

(s/^ . TN,- , , g^s to form an HII region. The MYSO phase is of particular 

^ , Massive stars (M> 8 Mo) form within dense molecular clouds of -^^^^^^^^ ^^^^^^^ ^^-^ ^^^^ ^^-^^ ^^^^ ^^^^ 

^ ' ^'l^f ' ^^1' Tl^^ '° [^P^ ^ 'u^i ''^ started core hydrogen burning, but major accretion is presum- 

■ ■ ■ observed direcdy before it has reached the main sequence. As a ^^^jj ■ ^ySOs are therefore a key stage for learning 

result, the formation mechanism and early evolution of inassive ^^^^^ ^^^^.^^ ^^^^^ ^^^^ a nd probing the various forma- 

stars IS still shrouded in mystery, with httle consensus in the field ^.^^ ^^^^^.^^ ^^^.^.^^^^ pj.^p^^gj (■Bernasconi & Maeder, 199^ 

as to a single formation theory, as opposed to that for low mass ^^^^^^ ^ ^^^^ 2003; Ballv & Zinnecker, 2005; Beuther et af 

star formation. |2qq7[) 

A tentative observationally based evolutionary sequence is ^, . ^^^o,^ u- u u i. j i 

... ■' . , , , , ■ ■ However, those MYSOs which have been widely 

emerging, starting with dense cores in the molecular cloud and , , ij- j j-.i j 

■ , , , , , , , . observed were mainly discovered serendipitously, and 

progressing through the hot molecular core and massive young ^, i. o -.i- m * o j . / 

.11 u- wiv/ivc/-,\ u iFu — i: — ril hnm k c . ii are often rather nearby. Sur ve ys w ith IRAS da ta (e.g . 

steUar object (MYSO) phases (lOur^wdll^ ICampb ell. Persson & MatthewsL & IPallaetall S 

the star begins to ionise the material around it and an ultra- „ ... ^^TTRTTTTTT T — j , • , ^ — c j 

iiT^ /iT^iiTT^ • i7 . • . . J ■ Sridhara n et aU 120021) have onlv considered bright, uncontused 

compact HII (UCHII) region forms, which then expands into an rn a o li. r. • i . j . j j 

. , , ,,TT • 1 .c 1 1 .o 11 1 IRAS sources. As such their samples tend to avoid the dense 

increasingly large HII region before the central star finally be- , r i i j . r ■ ■ i m a o 

. , . , -' parts of the galactic plane, due to confusion in the large IRAS 

comes visible in the optical as an OB star {; //ic-// t/ia// n u » • u »u • e 

... , ^ .jxrrT^ 1 1111 • 1 ■ J- 1 beam (45 x 240 at 12um), but this is where the maiority of 

We define an MYSO to be an embedded near or rmd-infrared u u i- a u • i j .vc j in-i^o/-. 

, . , . , , ^ T-. MYSOs should he. As such, the previously identified MYSOs 

(IR) source that IS luminous enough to be a young O or B-type i-i i » i. . i^.u i u i 

J J t- ^j-g unlikely to be representative of the class as a whole. 

* Based on observations collected at the European Southern The mid-infrared Galactic Plane su rvey with the Mid-course 

Observatory, La Silla, Chile (ESO Programmes 70.C-0069(A), 73.C- Space Experiment (MSX) by iPrice et al. (2001) has provided 

0314(A), 74.C-0267 and 077.C-0687(A)) better resolution data (-18" in all bands), allowing for most 

** E-mail:jcm@ast. leeds.ac.uk of these problems to be overcome. iLumsden et al.l (|2002|) de- 
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veloped colour selection criteria for MYSOs in the MSX point 
source catalogue (PSC), using the PSC fluxes. 2 Micron All- 
Sky Survey (2MASS) PSC fluxes (Strutskie et al., 2006) were 
used to help eliminate blue sources. These colour-cuts (Fg < 
Fi4 < F21, F2\IF% >2, Fa/F/c >5 and Fk/Fj >2), along with 
inspection to remove sources which are actually extended in 
the MSX observations, have returned a sample of ~2000 can- 
didate MYSOs. These are contaminated with both young and 
old dusty objects which have similar near and mid-IR colours 
(e.g. evolved stars, planetary nebulae (PN), HII regions). The 
Red MSX Source (RMS) Survey aims to remove the contam- 
inants from this MYSO candidate list via a series of ground 
based follow up ob servations at radio, mm and IR wavelengths 
dHoare et al.L l2005h . It will also provide data on the many new 
UCHII regions and evolved objects detected as a by-product of 
the survey process. 

At mid-IR wavelengths, the emission around MYSOs comes 
primarily from the warm (e.g. ~300K at lO/im) dust close to, 
and heated directly by, the young star. Observations of MYSOs 
at Ki2.2ijm) (e.g. iPorter et all Il998t iPreibisch et all l2003h 
show significant infrared excess from hot dust emission, sug- 
gesting that dust is present right up to the sublimation radius. The 
spherically symmetric mode ls of dust emission aroun d a 34M0, 
2.5xlO^L0 06 type star by IChurchwell et al.l (Il990h were ori- 
gionally intended for UCHII regions, but do not include Lyman 
a heating in an ionised zone and so are actually more ap- 
plicable to MYSOs. Their model 1 shows that the full-width 
at half maximum (fwhm) of intensity at 12yum (their Fig. 2a) 
is O.OOSpc (900AU), which at their distance of 2.6kpc corre- 
sponds to 0.36". This source lies at the near end of our dis- 
tance distribution and this fwhm is well below our derived limit 
for classification as a point source of 1.6" (see §3). The lu- 
minosity of the modelled source is also at the high end for 
our sample and since the emission size should scale approxi- 
mately as if it follows the size of Tj,,/,, the majority of our 
sources should be smaller than this. This basic situation will 
not change substantially in a disk-like geometry; most of the 
mid-IR emission will still come from the dense material close 
to the star, e.g. the models of Kessel, Yorke & Richling ( 1998) 
for a 8.4Mo star give a size of ~200AU at 12fj.m. However, the 
dust in the walls of the cavity can be heated to ~300K at fur- 
ther distances giving weak extended emission. While there are a 
few rare examples of resolved dust emission at ~lQfim around 
MYSOs (e.g. G35.2N De Buizer. 2006), the majority of obser- 
vations show only unres olved emission with the current gen- 
eratio n of facilities (e.g. iMarengo et all 120001 : iKraemer et al.L 
I2OOII; iDe Buizer etaTl l2005l) and we flierefore classify YSOs 
as point sources in the mid-IR. In the case of G35.2N we are 
probably viewing the system side on with respect to the disk, so 
the central source is too extincted to be visible. Were we to de- 
tect something similar we would expect the source to have quite 
a compact symmetric bipolar or monopolar appearance without 
strong radio continuum emission associ ated with it. 

For an UCHII region, the models of lHoare et al. I (|1991L their 
Fig. 7) suggest that 300K emission comes from the entire ionised 
region due to Lyman a heating. The mid-IR and radio continuum 
emission will therefore have similar a size and morphology for 
a given UCHII regions, and so will appear extended at 1" reso- 
lution at Ikpc, as will PN. Extended mid-IR objects can there- 
fore be removed from the MYSO candidate sample. While ra- 
dio continuum imaging can detect ionised regions, and may help 
identify unresolved UCHII regions or PN, it cannot detect the 
MYSOs which may lie near or superimposed upon HII regions. 
Proto-PN, low-mass stars between the post-AGB and PN phases. 



still have a low enough T^/f that they produce little radio contin- 
uum emission and can appear point-like at mid-IR wavelengths 
(Meixneret al., 1999). Evolved stars will also appear point-like, 
and so both will require identification using other information. 
In theory reddend RGB stars could be destinguished from YSOs 
using JHK colour-colour diagrams, but since sources selected 
by our colour selection criteria are either non-detections or have 
upper limit fluxes at J this is not possible in practise. In general, 
evolved stars usually appear isolated, rather than in clustered star 
formation regions and have little molecular CO, so both these 
facts can be used to diff'erentiate them from younger sources. 
Near-IR spectroscopy can also be used. 

High-resolution mid-IR imaging allows the combined flux 
of larger beam observations to be apportioned among the diff'er- 
ent sources if more than one object is present within the beam 
(e.g. MSX). This will improve the accuracy of any information 
derived from these data and allow the removal of MYSO can- 
didates from the list which are point-like but whose mid and 
far-IR fluxes in larger beam observations are actually dominated 
by extended emission from a nearby HII region. In addition, 
the improved astrometric information will be useful for subse- 
quent near-IR spectroscopy as the MSX coordinates are not good 
enough to identify the near-IR counterpart in crowded regions. 

Here we present mid-infrared observations undertaken in the 
southern hemisphere as part of the RMS survey. In §2 we report 
the details of the observational process. In §3 we present our re- 
sults, along with discussion of the astrometry (§3.1), photometry 
(§3.2), and a comparison with radio continuum data (§3.3). In §4 
we reach our conclusions. 



2. Observations 

Imaging observations were undertaken at IQAfim (Filter 
FWHM=1.02yL(m) with the Thermal Infrared Multi-mode 
Instrument 2 (TIMMI2; iKaufl et all |2003|) on the European 
Southern Observatory 3.6m telescope at La Silla, Chile over 4 
runs, totalling 12 nights between 2003 and 2006. The 10.4//m 
filter was chosen as it does not include either PAH emission fea- 
tures which could lead to extended emission, or the [Nell] line 
at 12.8/im which can be present in ionised regions. A nod-chop 
observing mode was used, with chopping perpendicular to nod- 
ding, a usual chop/nod throw of 15" and a chop frequency of 
6Hz. A pixel scale of 0.202"/pix was used, resulting in a field 
of view of 64"x 48" for each individual image. Standard stars 
were observed at ~2 hour intervals for flux calibration. For the 
majority of our observations, the seeing was no worse than 1.5" 
in the optical and close to the diffraction limit at lO/im. 

The blind pointing positional accuracy of the 3.6m telescope 
is 5", and so on the 10 nights of observations from 2004 onwards 
the following procedure was used to improve the positional ac- 
curacy of our observations. First an optically visible reference 
star, chosen to be <1° away from the target and bright enough 
to be detected quickly at 10/im (mv<12 and Fg >0.7Jy), was 
centred on a reference pixel. Then, the telescope was traversed 
by the known off'set between the reference star and the target co- 
ordinates. This ensured that the reference pixel in each object 
image was centred at the MSX target coordinates. The reference 
stars were chosen from the SAO star catalogue for the 2 nights of 
observations in 2004, but in some cases the reference stars were 
too faint. Therefore, a higher flux limit and the more extensive 
TYCHO catalogue was used for subsequent observations. This 
mechanism was not used for the 2 nights in 2003, and so we 
must assume that the targets observed in these observations are 
at the MSX positions. 
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A total of 346 targets were observed. The majority 
are outside the Spitzer Galactic Legacy Infrared Mid-Plane 
Survey Extraordinaire (GLIMPSE; Benjamin et al., 2003) re- 
gion (l=10°-65°, 295°-350°, |b|<r). While some of our obser- 
vations lie within this region, they were either observed before 
the GLIMPSE survey was released, or the GLIMPSE images 
are particularly difficult to interpret and so benefit from observa- 
tions at slightly better spatial resolution which are not affected by 
emission from transiently heated PAHs. Integration times were 
calculated in order to detect down to the faintest MSX fluxes 
(~0.1Jy) if all that flux was concentrated in a point source, with 
times ranging from 126s to 720s. 

Data reduction and analysis were performed using IDL 
scripts, with the reduction consisting of addition of the individ- 
ual image frames to form nod/chop pairs, addition of all the nod- 
chop pairs for one observation to form a total image and then 
mosaicing of the resultant images. 

3. Results 

Our results are presented as follows: first we report some general 
information on the observations as a whole. These are then fol- 
lowed by discussion of the table of data and figure of all targets 
available in the online version only. Data on the measured angu- 
lar size of the target objects are also presented. Next a more de- 
tailed discussion of the measurement and accuracy of the astro- 
metric (§3.1) and photometric (§3.2) data is provided. Finally, a 
comparison is made between our observations and current radio 
continuum data (§3.3) both from Walsh et al. ( 1998) and data ob - 
tained as part of the wider RMS efforts (Urquh art et~am2007ah . 

A total of 346 targets were observed, with the breakdown of 
detections and non-detections shown in Table 1 . With our chosen 
nod/chop throw of 15" we are able to detect sources up to this 
scale at least. A non-detection could therefore be due to the sur- 
face brightness of the object detected in the 18" MSX beam be- 
ing below our detection threshold. These objects are most likely 
either HII regions or PN, not MYSOs, and so are of less inter- 
est to our survey. Alternatively the flux at 9.7//m can be reduced 
by factors of 5 or more with respect to the 8//m continuum in 
parti cularly embedded YSOs due to the silicate absorption fea- 
ture jSmith et al.L 12000 ). This could lead to a non-detection, but 
should not be taken as proof of the absence of an MYSO. 

We differentiate between target regions showing single or 
multiple sources for several reasons. Firstly, evolved stars are 
unlikely to be near any other sources of mid-IR emission and 
so we would expect all targets showing multiple mid-IR emis- 
sion components to be related to star formation. Secondly, in 
cases where multiple components exist, the emission levels ob- 
served in larger beam measurements, such as IRAS and MSX, 
must be apportioned among these components. This may mean 
that extended emission near an MYSO candidate dominates, 
in which case the candidate may not be luminous enough to 
be considered an MYSO. Finally, this allows us to examine 
how many of our targets show potential for clustering and, in 
the case of point sources on or near extended emission, poten- 
tially triggered star formation. Early signs of star formation are 
often seen in the den se gas surrounding UCHII regions (e.g. 
iHofner & Churchwelll 1 19961) with good theoretical reasons for 
believing that their formation has been triggered by compression 
due to the expansion of the HII region (lElmegreeni ri998). 

The results of our astrometric, photometric and source size 
measurements for all targets are shown in Table 2 available in 
the online version only. We will further describe these measure- 
ments, and their general trends and implications below. Target 



Table 1. General Target Data. 





Object Type 


Number 


% of Total 


Detection 




303 


87.6 


Single 




218 


63.0 




Point (<L6") 


152 


43.9 




Extended (>1.6") 


66 


19.1 


Multiple 




85 


24.6 




Point (<1.6") 


41 


11.9 




Extended (>1.6") 


17 


4.9 




Extended & Point 


27 


7.8 


Non-Detections 


43 


12.4 



The criterion for distinguishing between point and 
extended sources is discussed in the §3. 



object names are given as the galactic coordinates of the MSX 
point source position. In the case of multiple sources, these 
are numbered with increasing distance from the MSX posi- 
tion. In the case where no object was detected in our obser- 
vations, a 3cr upper li mit flux h as been provided. The offset 
of the MSX PSC (V2.3 l Egan et a l., 2003) target position from 
the TIMMI2 position (column 5) and between the 2MASS PSC 
counterpart and the TIMMI2 coordinates for each object (col- 
umn 11), where a counterpart exists and is clear, are also pro- 
vided. A search for a potential 2MASS point source counterpart 
to each TIMMI2 source was conducted for angular separations 
of up to 5", though these are only likely to be true counterparts 
for angular separations of <2.4". The counterpart coordinates 
and offset from the TIMMI2 source are given in columns 9-11. 
Where astrometric information was not available, the primary 
object is assumed to be at the MSX coordinates. Nearby ra- 
dio continuum detections and non-detections, ei ther in our own 
interferometric observations with the ATCA dUrquhart et all 
2007a) at 3. 6cm and 6cm at a spa tial resolution of 1" or in 
the lit e rafiire (I Walsh et a l., 1998; Thor tipson. Urquhart & Whitel 
120041; IVan de Steene & Pottasch, 19931) . have been shown 
in column 12, with '-' indicating that no observation 
has been made for this target during either set of ob- 
servations. Ov e rall 31 8 of our tar gets we re observed by 
lUrquhart et all ( l2007ah and 1 1 by fWalsh et al. (1998)_, one 



(G336.4917-01.4741)bv lThompson. U rquhart & White (200 
and o ne (G344.9766-01.9071) by IVan de Steene & Pottasc 
(Il993h . leaving 15 with no radio continuum observation in ei- 
ther set. Of those sources with no radio observations, 12 are in 
the northern hemisphere where VLA observations have been ob- 
tained. In some cases, extended emission at 10.4yum does not 
have a radio continuum detection, probably due to the surface 
brightness of the radio continuum emission being below the de- 
tection limit and/or resolved out by the interferometry. 

Figure 1 (available online only) shows on the left the total 
field of view (79" x 63") of the mosaiced TIMMI2 images for 
all targets for which sources were detected. In cases w here radio 
emission has been detected bv Urquh art et al.l ilOOl d), their con- 
tours have been plotted on the TIMMI2 image, with preference 
given to the 6cm detections where both 6cm and 3.6cm obser- 
vations were made (e.g. G33 1.0309-1-01. 2056). The beam shape, 
size and position angle are shown in the lower left corner of the 
image. For targets with multiple (e.g. G260. 6877-01. 3930) or 
extended (e.g. G264. 291 8 -nO 1.4700) sources of emission, nega- 
tives and unphysical positives can lie near the real emission, due 
to the nodding and chopping required for mid-IR observations 
and subsequent mosaicing. These were removed by 'patching' 
the image using the STARLINK GAIA package PATCH tool to 
estimate the background levels from an annulus around the re- 
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gion to be patched. However, in regions with large amounts of 
extended emission, it was difficult to remove all the negatives, 
and in some cases difficult to tell what was real and what was 
not due to sources nodding or chopping onto other sources. In 
particular, it should be noted that while images and results for 
G291. 2725 -00. 7198 are presented here, they are tentative in that 
this region contains a large amount of extended 10.4yL(m emis- 
sion. Our aim is to show in this paper the full field of view with- 
out artifacts so that comparisons with other wavelengths, partic- 
ularly radio continuum maps, can be made easily. Any appar- 
ent structures other than those reported in Table 2 should not be 
over-interpreted. All the TIMMI2 data and images presented in 
this paper are included in the RMS survey databasqj- In partic- 
ular, the raw and mosaiced un-patched images are available, in 
addition to the patched mosaiced images, for those wishing to 
do a more in-depth comparison with their own data. The inten- 
sity scale of all images ranges from the median to the maximum 
level observed. For most images the scale is linear, but for those 
with both bright and faint sources on the same image a square 
root scale has been used. 

The right hand images in figure 1 show the K band 2MASS 
image for the target using a square-root intensity scale. The MSX 
target positions are indicated with a marking the 3cr posi- 
tional errors based on comparisons between the positions of our 
sources and the MSX coordinates for the whole sample (5.7", 
see §3.1). TIMMI2 detections are marked with an 'X' indicat- 
ing our 3cr astrometric accuracy (2.4"), derived from compar- 
isons of a subsample with associated 2MASS sources (again, see 
§3.1). For simple sources the TIMMI2 and MSX positions agree 
very well, but in more complex star formation regions, the MSX 
position can be 5-10" away from the dominant point source at 
10.4jum (e.g. G260.6877-01.3930). It should be noted that mid 
and near-IR sources do not always correspond to the same ob- 
ject, despite being close to each other since mid-IR objects can 
be heavily extincted in the near-IR. One should therefore not re- 
gard the nearness of sources of emission in TIMMI2 and 2MASS 
as an implicit indication that they come from the same source, 
though this may often be the case. 

Source sizes were measured using a 2-D gaussian fit to the 
source on the image, with the results for the major and minor 
axes shown in figure 2. In order to make a decision on whether an 
object is resolved or not, we examined the distribution of source 
sizes near the peak. The apparent size of an object is a combi- 
nation of object physical size, seeing and the diffraction limit of 
the 3.6m telescope (0.7" at lOyum). To remove the latter from 
our decision of what the maximum unresolved size of an object 
is, a gaussian fit was performed to the size distribution data af- 
ter a correction for the beam had been made using Osource-idlhf 

9 1 

^bean)^ ■ O^ly the data to the left of the peak was used in the fit, 
since all of these objects should be unresolved. The fit was then 
re-convolved to fit the measured source sizes. This returned a 
skewed gaussian with a slightly smaller spread on the lower side 
(cTi) than the higher side (0-2), since the addition in quadrature 
of Osource and 9beam-Q-T' has more effect on smaller values of 
Osource- Using 3cr2 above the mean from this fit to the major axis 
FWHM we define the maximum size of a point source as 1.6", 
as used in table 1. The difference between the peak of the major 
and minor FWHM distributions (1.23" and 1.07" respectively) 
is primarily due to the fact that the point spread function for 
TIMMI2 on the 3.6m telescope is slightly elliptical rather than 
spherical. This difference (0. 14") is actually less than the size of 
one pixel (0.202"). 

' The RMS database is accessible at |http://wvyw.ast.leeds.ac.uk/RMS| 



3.1. Astrometry 

Astrometry was carried out for all observations from 2004 on- 
wards as described in §2. The position of the centre of the object 
was obtained using a 2-dimensional gaussian fit to the source and 
the difference between this position and the reference pixel mea- 
sured. Some of the targets with multiple sources from the 2003 
observations were observed again in 2006 in order to provide ac- 
curate astrometry, which were then applied to the earlier obser- 
vations, which had longer exposure times. The offsets from the 
MSX target coordinates of all targets with astrometry are plotted 
in figure 3, with the offset for the primary source being used in 
cases where multiple sources are detected. The offsets are simi- 
lar between RA and Dec, and return an uncertainty (Icr) in the 
mean offset of 1 .9" in RA and 1 .7" in Dec (see Fig. 3). However, 
this assumes that the primary source in our observations corre- 
sponds to the MSX source, which is not always true. For exam- 
ple, in the case of the primary source for G34 1.1 28 1-00. 3466 
(see Fig. 1), MSX detects one slightly extended source while we 
see one bright and two fainter point sources near a region of ex- 
tended emission. The mean offset for both RA and Dec is not 
zero, but is small and well within the 2" quoted positional er- 
ror of the MSX catalogue. What is more, previous comparisons 
with the MSX PSC have also found a non-ze ro total offset (e.g. 
iLumsden et all 120021: lUrquhart et al.Ll2007al) . 

Let us now consider a subsample of the data for which we 
expect the TIMMI2 to MSX positions to be in good agreement. 
Sources which show only one unresolved source in our obser- 
vations and are free of extended emission in both the TIMMI2 
and MSX images are unlikely to have their positions skewed 
by diffuse emission or multiplicity. These sources are likely to 
be predominantly isolated evolved stars, which further reduces 
confusion. We find that the offsets for these objects are, on av- 
erage, slightly smaller than those for all targets (see Fig. 4.), 
returning an astrometric accuracy of 1.4" in RA and 1.1" in 
Dec. However, if we compare the TIMMI2 source positions with 
those for the source counterpart in the more accurate (0.2") 
2MASS PSC, where one exists, we get a much better correlation 
(0.8") than comparisons with MSX (see Fig. 4.). Since we have 
chosen this subsample in order to ensure that there is no confu- 
sion between counterparts in the different observations, we reach 
the conclusion that the offset between TIMMI2 and MSX posi- 
tions is dominated by the errors in the MSX positions. On the 
other hand, the offset between TIMMI2 and 2MASS coordinates 
gives a more accurate representation of the astrometric accuracy 
in our TIMMI2 observations (0.8"). The asymmetry seen in the 
offsets for this subsample between TIMMI2 and MSX are not 
seen when comparing between TIMMI2 and 2MASS, suggest- 
ing that this is due to systematics in the MSX positions and not 
those of our observations. 

3.2. Photometry 

Aperture photometry was performed on all detected objects. For 
isolated point sources, an aperture of 2.4" was used, with sky 
annulus inner and outer radii of 3" and 4" respectively, whilst 
for extended sources and those near other sources, these were set 
individually. In the case of a G326. 7249+00.6159 sources 1 & 2 
and G298. 1829-00.7860 sources 1 & 2, which have particularly 
small separations, aperture photometry was performed for the 
easiest source to separate and for both sources combined, with 
the flux for the second source obtained by subtracting the flux for 
the first source from the total flux. The errors in these fluxes were 
set based on the sky flux in the annulus of the smaller source. 
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Fig. 2. Major and minor FWHM size of detected targets. The vertical dashed-dot line in the left hand plot shows the 1 .6" boundary 
between point and extended sources. The dashed lines show skewed gaussian fits to the data. 
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Fig. 3. Offset of TIMMI2 astrometry from MSX coordinates, in both Ra and Dec, for all data with astrometry. The dashed lines 
show gaussian fits to the data. 



rather than the background levels in order to better reflect the 
uncertainty in these measurements. The standard observed near- 
est in time to a given object was used for calibration. 

The estimated error in the flux measurements was obtained 
by repeating the measurement procedure with the standard aper- 
ture in 4 different background regions and calculating the av- 
erage count levels. For frames with large amounts of extended 
emission, the error is likely to be larger as it is sometimes diffi- 
cult to find areas free from source emission. In the case of non- 
detections, we have included a 3cr upper limit based on the error 
calculated in this way. Comparison of standard star observations 
indicate that the overall photometric accuracy of our observa- 
tions is O.OSJy. 

Figure 5 shows the comparison between MSX 8jum (left) and 
12//m (right) fluxes and the derived TIMMI2 10. ^m flux for all 
detections, excluding the 18 targets which are MSX llfim non- 
detections. The distribution of points in the comparison with the 
MSX 12yum band appears slightly squashed at lower fluxes, but 
this is is due to the fact that this band has a sensitivity limit of 
~lJy, compared to the ~0.1Jy limit of the MSX Sjum band. The 
solid line is a least squares bisector fit to the data and the dashed 
lines indicate perfect correlation (i.e. y = x) in order to guide the 
eye. Therefore, the difference between the two lines shows the 
difference between the MSX and TIMMI2 fluxes. We do not ex- 
pect a perfect correlation between either MSX band fluxes and 



those obtained from our observations due to the larger beam 
size of the MSX observations. This results in extended emis- 
sion that is not within our typical 2.4" aperture being included 
within the MSX flux, while it is minimised in the TIMMI2 pho- 
tometry. Even for point sources we do not expect perfect agree- 
ment between MSX and TIMMI2 fluxes due to the differences 
in filter profiles and the rising slope of the spectral energy distri- 
bution towards longer wavelengths. The TIMMI2 filter includes 
the often deep 9.7jum absorption feature, and the MSX 8yum band 
contains some polycyclic-aromatic hydrocarbon (PAH) emission 
features, which GLIMPSE images (with far better spatial res- 
olution than MSX at a similar wavelength) of si milar regions 
show to be quite diffuse and extended (e.g. M17 iPovich et all 
2007*). Nevertheless, there is a reasonably good correspondence 
between the MSX and TIMMI2 fluxes, with the spread being 
larger at lower fluxes due to the fact that diffuse emission in the 
MSX images becomes more of an issue for fainter objects. 

In order to examine the effect of extended emission on the 
correlation between MSX and TIMMI2 fluxes, we again exam- 
ine the 'clean' subsample mentioned previously (§3.1), shown in 
figure 6. The solid line is a least squares bisector fit to the data 
and the dashed lines indicate perfect correlation (i.e. y = x) in 
order to guide the eye. From the better agreement between the 
fit to the data and the line of concordance, particularly for the 
comparison with the MSX 8/vm data, it is clear that the fluxes 
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Fig. 4. Offset of TIMMI2 astrometry from MSX (top) and 2MASS (bottom) coordinates, in both RA and Dec, for the 45 objects in 
the 'clean' subsample that have astrometric information and clear 2MASS PSC counterparts. The dashed lines show gaussian fits to 
the data. 




Fig. 5. Comparison of fluxes obtained from TIMMI2 observations with those from the MSX PSC for Syum (left) and 12/im (right) 
for all targets with 12fj.m detections. The solid line corresponds to a linear bisector least squares fit to the data in log-log space, 
while the dashed line is the line y = x in order to guide the eye. 



more closely agree for this subsample than for all observations. 
This is to be expected, since these sources are mostly evolved 
stars and other isolated sources with little diffuse emission. The 
remaining difference reflects the filter profile and spectral differ- 
ences already mentioned. The cut-off in the lower MSX 12/im 
flux is again due to the fact that the detector is less sensitive in 
this band than in the MSX S/im band. 



IWalsh et all (l200lh performed broad N (8-13jt/m) and Q 
(17.3-22.7jum) band observations of 8 of our targets (see table 
3). They observe higher fluxes (on average approximately twice) 
what we do, due to the broader N band filter chosen for their 
observations (as opposed to our 9.46-11.21//m filter) which in- 
cludes the PAH emission features at 8.7/im, more of the PAH 
feature at 11.3//m and is less dominated by the 9.7//m silicate 




Fig. 6. Comparison of fluxes obtained from TIMMI2 observations with those from the MSX PSC for the 57 targets in the 'clean' 
subsample with MSX llfim detections. The soUd Hne corresponds to a Hnear bisector least squares fit to the data in log-log space, 
while the dashed line is the line y = x in order to guide the eye. 



Table 3. Sources also observed bv lWalshetan(l2001h . 



Our Object 


Walsh etal. (2001) 


Name 


Source Number 


G259.9395-00.0419 


6,7 


G268.4222-00.8490 


11,12 


G269.4582-01.4713 


14 


G291. 2725-00.7198 


16 


G293. 8282-00.7445 


19 


G300.9674+01.1499 


21,22,23 


G3 13.4573+00. 1934 


36 


G345.0061+01.7944 


53 



absorption feature. Our astrometric positions agree with theirs 
to within their quoted error using guide stars of 7". 

3.3. Comparison with Radio Continuum Observations 

We compare our mid-IR imaging w ith the high resolution ( 1.5") 
radio continuum o bservations of lUrquhart et al.l (l2007al) and 
IWalsh et al] (Il998b . 83 (24%) have associated radio detections. 
In general, the agreement in alignment between our observa- 
tions and the radio data is quite good, with some interest- 
ing cases where resolved radio continuum emission seems to 
trace an extended HII region while the mid-IR emission reveals 
some additional unresolved sources (e.g. G298. 1829-00.7860 
and G329.421 1-00.1631). While one always has to be care- 
ful of projection effects, some of these may well be instances 
where the expanding HII region has triggered further star for- 
mation. The radio contours also show us that there are a few 
unresolved UCHII regions (e.g. G33 1.3546-1-01. 0638) and PN 
(e.g. G292.7329-00.2482) within our targets, but that these only 
account for ~7% and ~3% of all unresolved mid-IR sources re- 
spectively. This justifies the main presumption that point sources 
in star formation regions are primarily MYSOs. 

Within the ionised zone of a UCHII region, the gas is ionised 
by Lyman (A >912nm) continuum photons from the central star 
and emits thermal bremsstrahlung in the radio. Resonantly scat- 
tered Lyman a photons within the ionised zone provide the main 
source of heating for the dust within this region, which re-emits 
the energy at infrared wavelengths. The radio and infrared emis- 
sion should therefore be correlated. 



In order to investigate whether any correlation between ra- 
dio and mid-IR emission exists, we selected another subsample 
of sources. Starting from all sources where the radio and mid- 
IR sources seem to directly correspond, we eliminated objects 
which were substantially resolved (>8") in either the radio or 
mid-IR since the 6cm and 10.4/im fluxes are likely to be aff'ected 
by either being resolved out or due to nod/chop effects. We also 
eliminate sources which are optically thick at radio wavelengths 
from the spectral index between our 3cm and 6cm observations, 
as the radio flux will not reflect the total nebular gas emission. 
Due to the criteria that the object be optically thin, so that we 
are observing all of the emission in the radio, we also eliminated 
those objects where we do not have both 3cm and 6cm obser- 
vations. In addition, sources which are identified as PN within 
the RMS survey due to weak '^CO detections (Urq uhart et all 
|2007b), and/or their isolation and lack of extended emission in 
MSX and/or 2MASS images were eliminated from the subsam- 
ple. This resulted in a subsample of 23 objects. The 6cm and 
10.4yum fluxes for these sources are shown in figure 7, where 
we compare with fluxes tak en from the calculat ed spectra for 
dusty HII region models of [Arthur et all (|2004 their Fig. 8). 
Overall our results agree well with the range of fluxes expected 
for HII regions based on these m odels, and are also con sistent 
with the older simpler models of iNatta & Panagia (1976). The 
two points which do not seem to follow the observed trend are 
G268.6162-00.7389 (radio quiet) and G348.6972-01.0263 (ra- 
dio loud). 

G268. 6 162-00 .7389 is associated with a stellar group 
dPutra et all l2003h an d with the reflectio n nebula Bran 225. 
The near-IR spectra of 'Por ter et all ( Il998l) show Hel, Fell and 
OI emission lines in addition to hydrogen recombination lines 
but the line ratios are too low to be case B. The '^C O(J=1-0) 
distance measurements of I Wouterloot & Brandl ([1989) suggest 
that it is at a distance of 2kpc, resulting in Lf/jj=6682Lo from 
IRAS photometry. All this data supports the classification by 
Porter et al. ( 1998) of this source as a Herbig Be star, in which 
case the radio emission could be from a stellar wind. This is not 
in agreement with the optically thin spe ctral index of 0. 1 2+0.01 
returned from the radio observations of lUrquhartet alj (l2007al) 
which are more suggestive of a weak B type star HII region. 
However, the near-IR spectra are not consistent with that of a 
HII region. G348.6972-01.026 3 hes within part o f RCW122, 
and is definitely radio loud, since ! Walsh et al.l (Il998h detect even 
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Fig. 7. lOyum vs 6cm integrated radio flux for a subsample of sources. The dotted lin es correspond to the minimum and maximum 
Fecm/FiO/jjfi ratio from the calculated spectra for the dusty Hll models of [Arthur et alj f2004, their Fig. 8). The majority of points lie 
within or near the expected region. 



stronger radio emission than we do. ' Persson et al ] (ll976h calcu- 
late a silicate optical d epth at of Tg 7^,„ = 3.5 for this region and 
iMcBreen et all d 19851) calculate Ay=l 17 and T9.7^„,~6.5 for the 
whole cloud, so there is certainly significant absorption towards 
this source. This suggests that the silicate absorption feature, as 
well as general extinction may play a part in the lower than ex- 
pected 10.4/im flux of this source. 

4. Conclusions 

We have presented 10.4/im imaging observations for 346 candi- 
date MYSOs in the southern hemisphere, primarily outside the 
GLIMPSE survey region. We obtain photometric and astromet- 
ric accuracies of ~0.05Jy and 0.8" respectively. ~64% of our 
observations contain unresolved sources, which are most likely 
either MYSOs or evolved stars. 24% of our observations con- 
tain only sources of extended 10.4/im emission (FWHM>1.6"), 
which we expect primarily to be UCHll regions. For a subsam- 
ple of our sources with radio continuum data we examine the 
FecMi - Fi()^,„ relation, with our results being consistent with the- 
oretical predictions for UCHII regions. Our accurate astrometric 
data will aid near-lR spectroscopy in establishing object iden- 
tities, and make clear which near-lR source is the counterpart 
of that seen at mid and far-lR wavelengths. A small fraction of 
point sources (~7%) are unresolved UCHll regions, but these are 
readily identified by comparison with our radio continuum ob- 
servations. In cases where multiple sources of 10.4yL(m emission 
are detected (25%), our photometric data allow us to apportion 
the flux of larger beam observations between the different com- 
ponents, and so more accurately determine source luminosities. 

Mid-IR observations of all the northern RMS sources 
not within the GLIMPSE region have been completed using 
Michelle on UKIRT and Gemini-North, and will be the sub- 
ject of a future paper Preliminary identification of all RMS 
sources has been completed and is available on the database. 
Near-lR spectroscopic confirmation is ongoing, as is resolving 
the remaining near/far and multi-component '^CO kinematic 
distances. Work has also begun on obtaining far-lR luminosity 
information for RMS sources without good quality IRAS PSC 
fluxes. We are currently exploring using the IRIS and IGA data 
products as well as MIPSGAL where available. This will allow 
us to finally assign luminosities to all confirmed YSOs, allowing 
a final catalogue to be published. 
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